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Abstract 

The energy-loss formula of the production of gravitons by the binary is derived 
in the source theory formulation of gravity . Then, the quantum energy loss formula 
involving radiative corrections is derived. We postulate an idea that gravitational 
pulsars are present in our universe and that radiative corrections play a role in the 
physics of the cosmological scale. In the last part of the article, we consider so called 
electromagnetic pulsar which is formed by two particles with the opposite electrical 
charges moving in the constant magnetic field and generating the electromagnetic 
pulses. We think that the cosmological analogue is possible. 



1 PULSARS IN GENERAL 



Pulsars are the specific cosmological objects which radiates electromagnetic energy in the 
form of short pulses. They were discovered by Hewish et all. in 1967, published by Hewish 
et all. (1968) and specified later as neutron stars. Now, it is supposed that they are fast 
rotating neutron stars with approximately Solar mass and with strong magnetic fields 
(10 9 — 10 14 Gauss). The neutron stars are formed during the evolution of stars and they 
are the product of the reaction p + e~ — > n + v e , where the symbols in the last equation are 
as follows: proton, electron, neutron and the electron neutrino. The neutron stars were 
postulated by Landau (1932) at the year of the discovery of neutron in 1932. The pulsars 
composed from pions, or, hyperons, or, quarks probably also exist in universe, however it 
is not a measurement technique which rigorously determines these kinds of pulsars. 

Pulsars emit highly accurate periodic signals mostly in radio waves beamed in a cone of 
radiation centered around their magnetic axis. These signals define the period of rotation 
of the neutron star, which radiates as a light-house once per revolution. We know so called 
slow rotated pulsars, or, normal pulsars with period P > 20ms and so called millisecond 
pulsars with period P < 20ms. 

In 1974 a pulsar in a binary system was discovered by Hulse and Taylor and the 
discovery was published in 1975 (Hulse and Taylor, 1975). The period of rotation of 
normal pulsars increases with time and led to the rejection of the suggestion that the 
periodic signal could be due to the orbital period of binary stars. The orbital period of 
an isolated binary system decreases as it loses energy, whereas the period of a rotating 
body increases as it loses energy. 

The present literature concerns only the electromagnetic pulsars. The number of them 
is 1300 and they were catalogued with very precise measurements of their positions and 
rotation rates. The published pulse profiles are so called integrated profiles obtained by 
adding some hundred of thousands of individual pulses. The integration hides a large 
variation of size and shape from pulse to pulse of the individual pulsar. The radiation 
is emitted along the direction of the field lines, so that the observed duration of the 
integrated profile depends on the inclination of the dipole axis to the rotation axis, 
because it is supposed that the pulsar radiation is a radiation of the dipole in the magnetic 
field. It means that the radiation of these pulsars is the synchrotron radiation of charged 
particles. We know that the analogical ultrarelativistic charge moving in a constant 
magnetic field radiates the synchrotron radiation in a very narrow cone and such system 
can be considered as a free electron laser in case that the opening angle of the cone is very 
small. Of course the angle of emission of pulsars is not smaller than radians. In other 
words the observed pulsars are not the free electron lasers. The idea that pulsar can be 
a cosmical maser was also rejected. 

Pulsar radio emission is highly polarized, with linear and circular components. 
Individual pulses are often observed to be 100% polarized. The study of the polarization 
of pulsar is the starting point of the determination of their real structure. 

The only energy source of the pulsar is the rotational energy of the neutron star. The 
rate of the dissipation of the rotational energy can be determined. The moment of inertia 
is fairly accurately known from the theory of the internal structure and the rotational 
slowdown is very accurately measured for almost every pulsar. So most of energy is 
radiated as magnetic dipole radiation at the rotation frequency leading to a measure of 
the magnetic dipole and the surface field strength. 

The published articles on pulsars deal with the observation of the pulses and with 
the theoretical models. The observational results giving an insight into the behavior of 
matter in the presence of extreme gravitational and electromagnetic field are summarized 
for instance by Manchester (Manchester, 1992). The emission mechanism of photons are 



reviewed from a plasmatical viewpoint by Melrose (1992). The morfology of the radio 
pulsars is presented in the recent treatise by Seiradakis et. al. (2004). The review 
of properties of pulsars involving the radio propagation in the magnetosphere and of 
emission mechanism is summarized in the article by Graham-Smith (2003). At the same 
time there are, to our knowledge, no information on so called gravitational pulsars, or, on 
the models where the pulses are produced by the retrograde motion of the bodies moving 
around the central body. 

So, the question arises, if it is possible to define gravitational binary pulsar, where 
the gravitational energy is generated by the binary system, or, by the system where two 
components are in the retrograde motion. We suppose that in case of the massive binary 
system the energy is generated in a cone starting from the component of a binary and it 
can be seen only if the observer is present in the axis this cone. Then, the observer detects 
gravitational pulses when the detector is sufficiently sensitive. There are many methods 
for the detection of the gravitational waves. One method, based on the quantum states 
of the superfluid ring, was suggested by author (Pardy, 1989). 

We know that gravitational waves was indirectly confirmed by the observation of the 
period of the pulsar PSR 1913 + 16. The energy loss of this pulsar was calculated in 
the framework of the classical theory of gravitation. The quantum energy loss was given 
for instance by Manoukian (1990). His calculation was based on the so called Schwinger 
source theory where gravity is considered as a field theory of gravitons where graviton 
is a boson with spin 2, helicity ±2 and zero mass. It is an analogue of photon in the 
electromagnetic theory. 

In the following text we start with the source derivation of the power spectral formula 
of the gravitational radiation of a binary. Then, we calculate the quantum energy loss of 
a binary and the gravitational power spectrum involving radiative corrections. In the last 
part of an article, we consider so called electromagnetic pulsar which is formed by two 
particles with the opposite electrical charges which move in the constant magnetic field 
and generate the electromagnetic pulses. 

2 THE QUANTUM GRAVITY ENERGY LOOS OF 
A BINARY SYSTEM 

2.1 Introduction 

At the present time, the existence of gravitational waves is confirmed, thanks to the 
experimental proof of Taylor and Hulse who performed the systematic measurement of 
the motion of the binary with the pulsar PSR 1913+16. They found that the generalized 
energy-loss formula, which follows from the Einstein general theory of relativity, is in 
accordance with their measurement. 

This success was conditioned by the fact that the binary with the pulsar PSR 1913+16 
as a gigantic system of two neutron stars, emits sufficient gravitational radiation to 
influence the orbital motion of the binary at the observable scale. 

Taylor and Hulse, working at the Arecibo radiotelescope, discovered the radiopulsar 
PSR 1913+16 in a binary, in 1974, and this is now considered as the best general 
relativistic laboratory (Taylor, 1993). 

Pulsar PSR 1913+16 is the massive body of the binary system where each of the 
rotating pairs is 1.4 times the mass of the Sun. These neutron stars rotate around each 
other with a period 7.8 hours, in an orbit not much larger than the Sun's diameter. Every 
59 ms, the pulsar emits a short signal that is so clear that the arrival time of a 5-min 
string of a set of such signals can be resolved within 15 /xs. 



A pulsar model based on strongly magnetized, rapidly spinning neutron stars was 
soon established as consistent with most of the known facts (Huguenin et al, 1968); 
its electrodynamical properties were studied theoretically (Gold, 1968) and shown to 
be plausibly capable of generating broadband radio noise detectable over interstellar 
distances. The binary pulsar PSR 1913+16 is now recognized as the harbinger of a 
new class of unusually short-period pulsars, with numerous important applications. 

Because the velocities and gravitational energies in a high-mass binary pulsar system 
can be significantly relativistic, strong-field and radiative effects come into play. The 
binary pulsar PSR 1913+16 provides significant tests of gravitation beyond the weak- 
field, slow-motion limit (Goldreichet al., 1969; Damour et al., 1992). 

The goal of this section is not to repeat the derivation of the Einstein quadrupole 
formula, because this has been performed many times in general relativity and also in the 
Schwinger source theory in the weak- field limit (Manoukian, 1990). 

We show that just in the framework of the source theory it is easy to determine 
the quantum energy-loss formula of the binary system. The energy-loss formula can be 
generalized in such a way it involves also the radiative corrections. 

Since the measurement of the motion of the binaries goes on, we hope that future 
experiments will verify the quantum version of the energy-loss formula, involving also the 
radiative corrections. 

2.2 The source theory formulation of the problem 

We show how the total quantum loss of energy caused by the production of gravitons, 
emitted by the binary system of two point masses moving around each other under their 
gravitational interaction, can be calculated in the framework of the source theory of 
gravity. 

Source theory (Schwinger, 1970, 1973, 1976) was initially constructed to describe the 
particle physics situations occurring in high-energy physics experiments. However, it was 
found that the original formulation simplifies the calculations in the electrodynamics and 
gravity, where the interactions are mediated by photon and graviton respectively. The 
source theory of gravity forms the analogue of quantum electrodynamics because, while in 
QED the interaction is mediated by the photon, the gravitational interaction is mediated 
by the graviton (Schwinger, 1976). 

The source theory of gravity invented by Schwinger is linear theory. So, the question 
arises if it is in the coincidence with the Einstein gravity equations which are substantially 
nonlinear. The answer is affirmative, because the coincidence is only with the linear 
approximation of the Einstein theory. The experimental results of the Schwinger theory 
are also in harmony with experiment. The quadrupole formula of Einstein also follows 
from the Schwinger version. 

The unification of gravity and electromagnetism is possible only in the Schwinger 
source theory. It is possible when, and only when it is possible the unification of forces. 
And this is performed in the Schwinger source theory of all interactions where force is 
of the Yukawa form. The problem of unification is not new. We know from the history 
of physics that the Ptolemy system could not be unified with the Galileo-Newton system 
because in the Ptolemy system it is not defined the force which is the fundamental quantity 
in the GN system and the primary cause of all phenomena in this system. 

Einstein gravity uses the Riemann space-time where the gravity force has not the 
Yukawa dynamical form. The curvature of space-time is defined as the origin of all phe- 
nomena. The gravity force in the Einstein theory is in the antagonistic contradistinction 
with the Yukawa force in the quantum field theory and therefore it seems that QFT, 
QED, QCD and EL.- WEAK theory cannot be unified with the Einstein gravity. 



Manoukian (1990) derivation of the Einstein quadrupole formula in the framework of 
the Schwinger source theory is possible because of the coincidence of the source theory 
with the linear limit of the Einstein theory. 

Our approach is different from Manoukian method because we derive the power 
spectral spectral formula P{oj) of emitted gravitons with frequency uj and then using 
relation (—dE/dt) = J duP(u) we determine the energy loss E. In case of the radiative 
correction, we derive only the power spectral formula in the general form. 

The mathematical structure of P(w) follows directly from the action W and while 
in the case of the gravitational radiation, the formula is composed from the tensor of 
energy-momentum, then, in the case of the electromagnetic radiation formula, it involves 
charged vector currents. 

The basic formula in the source theory is the vacuum-to- vacuum amplitude (Schwinger 
et al., 1976): 

(0 + |0_> = e^ 5 \ (1) 

where the minus and plus tags on the vacuum symbol are causal labels, referring 
to any time before and after region of space-time, where sources are manipulated. 
The exponential form is introduced with regard to the existence of the physically 
independent experimental arrangements, which has the simple consequence that the 
associated probability amplitudes multiply and the corresponding W expressions add 
(Schwinger et al., 1976; Dittrich, 1978). 

In the flat space-time, the field of gravitons is described by the amplitude (1) with the 
action (Schwinger, 1970) (c = 1 in the following text) 

W(T) =4nG J{dx){dx') 

(2) 



x 



T^(x)D+(x - x')T^(x') - -T(x)D + (x - x')T(x' 



1, 
2" 

where the dimensionality of W(T) is the same as the dimensionality of the Planck constant 
h; T^y is the tensor of momentum and energy. For a particle moving along the trajectory 
x = x(i), it is defined by the equation (Weinberg, 1972): 

T^(s) = ^(x-x(f)), (3) 
where is the relativistic four-momentum of a particle with a rest mass m and 

p" = (£,p) (4) 

p% = -m 2 , (5) 
and the relativistic energy is defined by the known relation 

where v is the three- velocity of the moving particle. 

Symbol T(x) in formula (2) is defined as T = g^T^, and D + (x — x') is the graviton 
propagator whose explicit form will be determined later. 

The action W is not arbitrary because it must involve the attractive force between the 
gravity masses while in case of the electromagnetic situation the action must involve the 
repulsive force between charges of the same sign. It is very surprising that such form of 
Lagrangians follows from the quantum definition of the vacuum to vacuum amplitude. It 



was shown by Schwinger that Einstein gravity also follows from the source theory, however 
the method of derivation is not the integral part of the source theory because the source 
theory is linear and it is not clear how to establish the equivalence between linear and 
nonlinear theory String theory tries to solve the problem of the unification of all forces, 
however, this theory is, at the present time, not predictable and works with so called 
extra-dimensions which was not observed. It is not clear from the viewpoint of physics, 
what the dimension is. It seems that many problems can be solved in the framework of 
the source theory. 

2.3 The power spectral formula in general 

It may be easy to show that the probability of the persistence of vacuum is given by the 
following formula (Schwinger et al., 1976): 

|(0+|0_)| 2 = expj--lmwj = expj-y dtdcu^-P(uj, t) J , (7) 

where the so-called power spectral function P(uj,t) has been introduced (Schwinger et 
al., 1976). In order to extract this spectral function from Im W, it is necessary to know 
the explicit form of the graviton propagator D + (x — x'). The physical content of this 
propagator is analogous to the content of the photon propagator. It involves the gravitons 
property of spreading with velocity c. It means that its explicit form is just the same as 
that of the photon propagator. With regard to the source theory (Schwinger et al., 1976) 
the ^-representation of D + (x) in eq. (2) is as follows: 



D + (x-x>)=J^e ik ^D(k), 



where 



which gives 



D + (x -x') = — dcu—^ —L e -">\t-t I. 10 

47T 2 JO X-X' 



Now, using formulas (2), (7) and (10), we get the power spectral formula in the 
following form: 



P(u,t) =4nGu J (ebz)(dx')dt! Sm( ^* * ! cosu{t - t' 



x 



T^(x,f)T^(x',0 - i^T^(x,0^T^(x',t') 



2.4 The power spectral formula for the binary system 



In the case of the binary system with masses mi and ttt.2, we suppose that they move in a 
uniform circular motion around their center of gravity in the xy plane, with corresponding 
kinematical coordinates: 

xi(t) = ri(icos(cj £) +jsin(cj £)) (12) 



X2(0 = r 2 (i cos(c<;ot + n) + j sin(cj £ + 7r)) 



(13) 



with 

Vi(t) = dxi/dt, uj = Vi/ri, v { = |v;| (i = 1, 2). (14) 
For the tensor of energy and momentum of the binary we have: 

T^(x) = ^*(x - Xl (f)) + ^5(x - x 2 (f)), (15) 

Cj\ Hi 2 

where we have omitted the tensor t^ u , which is associated with the massless, gravitational 
field distributed all over space and proportional to the gravitational constant G (Cho et 
al., 1976): 

After insertion of eq. (15) into eq. (11), we get: 

Ptotal(w,t)=Pl(w,t)+P 12 (w,t)+P 2 (uJ,t), (16) 

where (t' — t — r): 

Gu sin[2cjri sin(cu r/2)] 

Pi{u,t) = / dr — — cos^r 

ri7TJ-oo sin (ujftT/2) 



x ^(wgrjcoswor-l) 2 -^, (17) 

Gcv [°° sin[2^r 2 sin(wor/2)] 

P 2 (w,f) = / dr — — cos^r 

r 2 7r J-oo sm(a;oT/2) 

E 2 2 K 2 r 2 2 cos^ r - l) 2 - , (18) 

„ , . AGuj r°° , sin w[r 2 + ra + 2r x r 2 cos(w T")l 1 ^ 2 
P 12 iuj,t) = / dr — ^ 5 — —-jz — cos lot 

v ' it J-oo [r\ + r\ + 2rir 2 cos(cj r)] 1 / 2 

x (^(c^cos^r + l) 2 - . (19) 

2.5 The quantum energy loss of the binary 

Using the following relations 

lo t = tp + 2nl, ¥>e(-7r,7r), / = 0, ±1, ±2, ... (20) 

i=oo oo 

cos27r/ — = cj 5(^ — cu /), (21) 

l=— oo ^0 i=— oo 

we get for Pi(to,t), with a; being restricted to positive: 

oo 

P(cu, = £ % " "o0^("> 0- (22) 
i=i 

Using the definition of the Bessel function J 2 i(z) 

■hi^z) — — I dip cos fzsin^ cos/y?, (23) 
27T j-tt V 2 / 

from which the derivatives and their integrals follow, we get for Pu and P 2 i the following 
formulas: 



4 



2Guu 



(im-v^-^) i dxj 2l { X ) 



+ 4£ 2 (1 - v^vjj'^vd) + AE%vUZ$Vil)), 



1, 2. 



(24) 



Using r 2 = r\ + e, where e is supposed to be small in comparison with radii r\ and r 2 , 
we obtain 



[r 2 + r 2 + 2rir 2 cos y?] 1 / 2 ~ 2a cos , 



with 



a = r l( l + ->. 



(25) 



(26) 



So, instead of eq. (19) we get: 



IGoj f°° sin[2a;acos(a;or/2)] 



P 12 (uj,t) = / (I, : --T- — COS a,' 7" 

air J-oo cos(a;or)/ 



cos(cj t)/2] 
x I E 1 E 2 (u 2 r 1 r 2 cos lu t + l) 2 



(27) 

\ 2E\E 2 J 

Now, we can approach the evaluation of the energy-loss formula for the binary from 
the power spectral formulas (24) and (27). The energy loss is defined by the relation 

dU 



au f ^, x , 



J du^6(u - wol)Pu + / Pi2Hdiu = -f f (Ui + u 2 + U 12 ). 



Or, 



I i,l 

From Sokolov and Ternov (1983) we learn Kapteyn's formulas: 



(28) 



(29) 



J22U' 2l (2lv) = — 
i=i \ l 



2\2 ' 



V 2 ) 



and 



00 r2lv 



v 2 f 



The formula M 2 [(2lv) = can be obtained from formula 

1 7/ ^ , 1 



(30) 



(31) 



(32) 



by its differentiation with the respect to v (Schott, 1912). 

Then, after application of eqs. (30), (31) and (32) to eqs. (24) and (28), we get: 



dUi 2Gtu 



dt 



Hh- \ v l 2E 2 v 3 + AE 2 v 4 

2E 2 3(1 -^Y 11 11 



(33) 



Instead of using Kapteyn's formulas for the interference term, we will perform a direct 
evaluation of the energy loss of the interference term by the ^-integration in (27). So, 
after some elementary modification in the u;-integral, we get: 



dU 



12 



dt 



/oo rco 
dr / duuje~ lU}T sin[2u;a cos ojqt] 
-co J —CO 



rco 

= / P{u)duj = 
Jo 

B(C cos ujqt + 1) 2 - D 



cos(cj t/2) 



with 



G 



A = —, B = E 1 E 2 , C = Vl v 2 , D = !!^A, 
avr' 2 ' 2 ' 2E 1 E 2 



(34) 



(35) 



Using the definition of the 5-function and its derivative, we have, instead of eq. (34), 
with v = aujQ. 



dU 



12 



dt 



°° [B(Ccosx + l) 2 - D] 

(tCC ~ ■ r X 

oo cos(x/2) 



[S'{x - 2v cos(x/2)) - 5'{x + 2v cos(x/2))] 



Putting 



x — 2v cos(x/2) = t, 



in the first <5'-term and 



y + 2vcos(y/2)=t, 
in the second 5'-term, we get eq. (36) in the following form: 

dU 12 



(36) 



(37) 



(38) 



dt 



/CO 
dtS'(t) x 



f [B(C cosx + l) 2 - D] [BjCcosy + l) 2 - D] ' 
\(x-t)(l+vsm(x/2)) ~ (y + t)(l-vsm(y/2))\ 

Using the known relation for a 5-function: 

J dtf(t)5'(t) = -f (0), 
we get the energy loss formula for the synergic term in the form: 



(39) 



(40) 



dU 12 d j [B(C cosx + l) 2 — D] [B(Ccosy+ l) 2 -D] \ 

dt U ° V71 di\(x -t)(l + vsm(x/2)) ~ (y-t)(l -vsm(y/2)) J 

and we recommend the final calculation of the last formula to the mathematical students. 

Let us remark finally that the formulas derived for the energy loss of the binary (33) 
and (41) describes only the binary system and therefore their sum has not the form of 
the Einstein quadrupole formula. The sum forms the total produced gravitational energy, 
and involves not only the radiation of the individual bodies of the binary, but also the 
interference term. The problem of the coincidence with the Einstein quadrupole formula 
is open. 



(41) 



3 THE POWER SPECTAL FORMULA IVOLVING 
RADIATIVE CORRECTIONS 



3.1 Introduction 

We here calculate the total quantum loss of energy caused by production of gravitons 
emitted by the binary system in the framework of the source theory of gravity for the 
situation with the gravitational propagator involving radiative corrections. 

We know from QED that photon can exist in the virtual state as the two body system 
in the form of the electron positron pair. It means that the photon propagator involves 
the additional process: 

7 -> e + + e~ -> 7. (42) 

In case of the graviton radiation, the situation is analogical with the situation in QED. 
Instead of eq. (42) we write 

g _> 2e + + 2e" -> g, (43) 

where g is graviton and number 2 is there in order to conserve spin also during the virtual 
process. 

Equation (43) can be of course expressed in more detail: 

g -> 7 + 7 -»■ (e + + O + (e + + e") -»• 7 + 7 -»■ 5- (44) 

We will show that in the framework of the source theory it is easy to determine 
the quantum energy loss formula of the binary system both in case with the graviton 
propagator with radiative corrections. 

We will investigate how the spectrum of the gravitational radiation is modified if we 
involve radiation corrections corresponding to the virtual pair production and annihilation 
in the graviton propagator. Our calculation is an analogue of the photon propagator 
with radiative corrections for production of photons by the Cerenkov mechanism (Pardy, 
1994c,d). 

Because the measurement of motion of the binaries goes on, we hope that the future 
experiments will verify the quantum version of the energy loss formula following from the 
source theory and that sooner or later the confirmation of this formula will be established. 



3.2 The binary power spectrum with radiative corrections 

According to source theory (Schwinger, 1973; Dittrich, 1978; Pardy, 1994c, d), the photon 
propagator in the Minkowski space-time with radiative correction is in the momentum 
representation of the form: 

D(k)=D(k) + 5D(k), (45) 



D(k) ' 



k| 2 - (A; ) 2 - ie 
a(M 2 ' 



km 2 |k| 2 -(fc°) 2 + K£-ie y ! 



where m is mass of electron and the last term in equation (44) is derived on the virtual 
photon condition 

\k\ 2 -(k°) 2 = -^-. (47) 

The weight function a(M 2 ) has been derived in the following form (Schwinger, 1973; 
Dittrich, 1976): 

a 1 f, 2m2 \ f, 4m 2 \ V2 
< M) = T.W\ 1+ Mt){ l -Mt) • (48) 

We suppose that the graviton propagator with the radiative corrections forms the 
analogue of the photon propagator. 

Now, with regard to the definition of the Fourier transform 

D + (x- X ') = J^e^-^D(k), (49) 

we get for 5D + the following relation (c = H — 1): 

j roo 

5D + (x -x') = — dM 2 a(M 2 ) 

AlT 2 Jim 2 



X / duj ^ : LL e -^l*-*'l. (50) 

J |x — x'| 

The function (50) differs from the gravitational function " D + n in (9) especially by the 
factor 

[u 2 -M 2 ) 1 ' 2 (51) 

in the function 'sin' and by the additional mass-integral which involves the radiative 
corrections to the original power spectrum formula. 

In order to determine the additional spectral function of produced gravitons, corre- 
sponding to the radiative corrections, we insert D + {x — x') + 5D + (x — x') into eq. (2), 
and using eq. (11) we obtain (factor 2 from the two photons is involved): 

6P(u,t) = — / (dx)(dx')dt f / dM 2 a(M 2 ) 

7T J J Am 2 

sinfL^-MWx-x'l} 

x - , -7T ^cosuit-t') 

|x — x'| 

x T^(x, t)^ Q ^T^(x', if) - ^T^x, t) 9afi T^(x!, If)] . (52) 
Then using eqs. (16), (17), (18) and (19), we get 

8P total {u, t) = SP^u, t) + 5P 2 (lu, t) + 5P 12 (u, t), (53) 

where (f — t — t): 

rrW x 2Glu r°° , f°° „, 2 / _ 2x sin{2(cj 2 - M 2 ) 1/2 risin(cj r/2)} 

SP 1 (u,t) = / dr dM 2 a(M 2 ) — ^ / -\ — v 1 n coscur 

r\-K J-oo Jim 2 sin(o;or/2) 



x ^KMcoscor-l) 2 -^, (54) 

rrW x 20j /-oo i-oo / _ 2x sin{2(cj 2 - M 2 ) 1/2 r 2 sin(cj r/2)} 

5P 2 (^,*) = / dr dM 2 a(M 2 ) — ^ / — v 1 Jf cosuT 

r 2 ir J-oo Jim? sm(u T/2) 

x ^(^cos^or-l) 2 -^, (55) 



2E% 



rsi -ill) r°° f°° 

5P 12 {uj,t) = — dr dM 2 a{M 2 

7T J-oo JAm? 

sin{(g> 2 - M 2 ) 1/2 [r 2 + rf + 2r 1 r 2 cos(iuqt)] 1 / 2 } 
[r\ + r\ + 2rir 2 cos(cjqt)] 1 ^ 2 



cos cur 



x ^i^K^a cosc^t + l) 2 - ^r^rj • (56) 

The explicit determination of the power spectrum is the problem which was solved 
by author in 1994. The solution was performed only approximately. Here also can be 
expected only the approximative solution. From this solution can be then derived the 
energy loss as in the previous article. 

Let us show the possible way of the determination of the spectral formula. If we 
introduce the new variable s by the relation 

uj 2 — M 2 = s 2 ; -dM 2 = 2sds (57) 
then, instead of eqs. (54), (55) and (56) we have 

. 2Guj f°° , r s * . , . . 9 n. sin{2sri sin(u;o^/2)i 

8PAu,t) = / dr / (2sds)a(u 2 - s 2 ) i v " 1 js coslut 

V ' ' J-oo J S1 V ' K ' sin cunT 2) 



Tin J-oo J si sin(c<JoT/2) 

2E\ 

2Gu f°° r s 2 2 2 sin{2sr 2 sin(a;or/2)} 



x E 2 K 2 r 2 cosc r-l) 2 -^ , (58) 



<JP 2 (w,0 = / dr / 2sds Wu; 2 - s 2 — 1 — 

V ' 7 r 2 vr J-oo ; V ; sin(c; r/2) 



COSCJT 



x l^w^coswor-l)'-^!, (59) 



8G<j yoo r-2 sin{2s[r 2 + r 2 + 2nr 2 cos(uqt)] 1 / 2 } 

dPi 2 (cj,t) = / dr / (2sds)a(ur - s z ) 5 — — — - 

7T 7-oo 7 Sl v [r 2 + r| + 2rir 2 cos(cj r)] 1 / 2 



2 2 



where 



x cos^r I E x E 2 {^r x r 2 cosu; t + l) 2 - 1 , (60) 



si = c<j 2 — 4m 2 , s 2 = oo. (61) 



It seems that the rigorous procedure is the r-integration as the first step and then 
s-integration as the second step (Pardy, 1994c, d). While, in case of the linear motion 



the mathematical operations are easy (Pardy, 1994c), in case of the circular motion there 
are some difficulties (Pardy, 1994d). The final form of eq. (60) is recommended for the 
mathematical experts. 

The energy loss is given as follows: 



Now, let us go to the discussion on the electromagnetic system of the two opposite 
charges moving in the constant magnetic field and producing the pulse synchrotron 
radiation. 



4.1 Introduction 

Here, the power spectrum formula of the synchrotron radiation generated by the electron 
and positron moving at the opposite angular velocities in homogeneous magnetic field is 
derived in the Schwinger version of quantum field theory. 

It is surprising that the spectrum depends periodically on radiation frequency u and 
time which means that the system composed from electron, positron and magnetic field 
forms the pulsar. 

We will show that the large hadron collider (LHC) which is at present time under 
construction in CERN can be considered in near future also as the largest electromagnetic 
terrestrial pulsar. We know that while the Fermilab's Tevatron handles counter-rotating 
protons and antiprotons in a single beam channel, LHC will operate with proton and 
proton beams in such a way that the collision center of mass energy will be 14 TeV and 
luminosity 10 34 cm _1 s -2 . To achieve such large luminosity it must operate with more than 
2800 bunches per beam and a very high density of particles in bunches. The LHC will 
also operate for heavy Pb ion physics at a luminosity of 10 27 cm~ 1 s~ 2 (Evans, 1999). 

The collision of particles is caused by the opposite directional motion of bunches. Or, 
if one bunch has the angular velocity u, then the bunch with antiparticles has angular 
velocity —uo. Here we will determine the spectral density of emitted photons in the 
simplified case where one electron and one positron move in the opposite direction on a 
circle. We will show that the synergic spectrum depends periodically on time. This means 
that the behavior of the system is similar to the behavior of electromagnetic pulsar. The 
derived spectral formula describes the spectrum of photons generated by the Fermilab 
Tevatron. In case that the particles in bunches are of the same charge as in LHC, then, 
it is necessary to replace the function sine by cosine in the final spectral formula. Now 
let us approach the theory and explicit calculation of the spectrum. 

This process is the generalization of the one-charge synergic synchrotron-Cerenkov 
radiation which has been calculated in source theory two decades ago by Schwinger et al. 
(1976). We will follow the Schwinger article and also the author articles (Pardy, 1994d, 
2000, 2002) as the starting point. Although our final problem is the radiation of the two- 
charge system in vacuum, we consider, first in general, the presence of dielectric medium, 
which is represented by the phenomenological index of refraction n and it is well known 
that this phenomenological constant depends on the external magnetic field. Introducing 
the phenomenological constant enables to consider also the Cerenkovian processes. Later 
we put n = 1. 

We will investigate here how the original Schwinger (et al.) spectral formula of the 
synergic c synchrotron Cerenkov radiation of the charged particle is modified if we consider 
the electron and positron moving at the opposite angular velocities. This problem is an 




(62) 
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analogue of the linear (Pardy, 1997) and circular problem solved by author (Pardy, 2000). 
We will show that the original spectral formula of the synergic synchrotron-Cerenkov 
radiation is modulated by function 4sin 2 (wt) where uj is the frequency of the synergistic 
radiation produced by the system and it does not depend on the orbital angular frequency 
of electron or positron. We will use here the fundamental ingredients of Schwinger source 
theory (Schwinger, 1970, 1973; Dittrich, 1978; Pardy, 1994c, d, e) to determine the power 
spectral formula. 



4.2 Formulation of the electromagnetic problem 

The basic formula of the Schwinger source theory is the so called vacuum to vacuum 
amplitude: (0+ 1 0_) = exp{|VF}, where in case of the electromagnetic field in the medium, 
the action W is given by the following formula: 

W = 7^J (dx)(dx')JV(x)D +llu (x - x')r(x'), (63) 

where 

D% = ^[ g ^ + (1 - n' 2 )^(3 v ]D + (x - x'), (64) 

where (3^ = (1, 0), J M = (eg, J) is the conserved current, /x is the magnetic permeability 
of the medium, e is the dielectric constant od the medium and n = ^JTjl is the index of 
refraction of the medium. Function D + is defined as in eq. (10) (Schwinger et al., 1976): 

D + (x - x') = -V [°° d /"^ 1 *;*^ !*-*'! (65) 
4n 2 cJo |x — x'| 

The probability of the persistence of vacuum follows from the vacuum amplitude (1) 
where Im W is the basis for the following definition of the spectral function P(u,t): 

_2 Im w A _ f dtduJ E^A. ( 66 ) 
n J nuo 

Now, if we insert eq. (64) into eq. (66), we get after extracting P(u,t) the following 
general expression for this spectral function: 



P(u,t) = / dxdx'dt' 

An 2 n 2 J 



sin^lx- x' 

cj 



n 2 



X COs[cj(t-t / )][^(x,t)^(x / ,t / ) - -rJ(x,t) • J(x',t')], ( 67 ) 

c 2 

which is an analogue of the formula (11). 

Let us recall that the last formula can be derived also in the classical electrodynamic 
context as it is shown for instance in the Schwinger article (Schwinger, 1949). The 
derivation of the power spectral formula from the vacuum amplitude is more simple. 



4.3 The radiation of two opposite charges 

Now, we will apply the formula (67) to the two-body system with the opposite charges 
moving at the opposite angular velocities in order to get in general synergic synchrotron- 
Cerenkov radiation of electron and positron moving in a uniform magnetic field 

While the synchrotron radiation is generated in a vacuum, the synergic synchrotron- 
Cerenkov radiation can produced only in a medium with dielectric constant n. We suppose 



the circular motion with velocity v in the plane perpendicular to the direction of the 
constant magnetic field H (chosen to be in the +z direction). 

We can write the following formulas for the charge density g and for the current density 
J of the two-body system with opposite charges and opposite angular velocities: 



and 



with 



£>(x, t) = e5(x. — x.i(t)) — e5(x — x 2 (t)) 



J(x,t) = evi(f)<$(x-xi(i)) -ev 2 (f)<J(x-x 2 (i)) 



(68) 



(69) 



xi(£) = x(i) = R(i cos(cjq^) + j sin(u;ot)), 



(70) 



x 2 (£) = R(icos(—ujot) + jsin(— uot) = x(— uo,t) = x(— t). 



(71) 



The absolute values of velocities of both particles are the same, or |vi(t)| = |v 2 (i)| = v, 
where (H = |H|, E = energy of a particle) 



BE 

v(t) — dx/dt, iv — v/R, R = ——, B = v/c, v — Ivl 

en 



(72) 



After insertion of eqs. (68)-(71) into eq. (67), and after some mathematical operations 
we get 

P(u,t) = - — ^J^ dt'cos(t-t>) £ ("ir 



X 



1 _ VjgWg) w2 



*J=1 

sin^lx^-x^o 



(73) 



Let us remark, that for situation of the identical charges, the factor (— must be 
replaced by 1. 

Using t' = t + r, we get for 



Xj(t) x j(t ) — A^, 



| | = [R 2 + R 2 - 2RRcos(lu t + a i:j )] 1/2 = 2R 
where were evaluated as follows: 



sin 



ujqt + a 



'■j 



an = 0, «i2 = 2u t, a 2 i = 2u t, a 2 2 = 0. 

Using 

Vi(t) ■ Vj(t + r) = uj 2 R 2 cos(uj t + aiij), 
and relation (75) we get with v = u R 

2 



(74) 
(75) 

(76) 
(77) 



P(u,t) = 



LU /X 2 



An 2 n 2 



e 2 dr cos our 



x 



1 COs(iO T + Otij) 



sin 



2Rsm( (uJ ° T + aij) ) 



1. 



(78) 



Introducing new variable T by relation 



for 



U T + Oiij = u T 

every integral in eq. (78), 



we get P{uj,t) in the following form 

-^^^ r d T t 



_ -l) i+J 



x cos(cjT CKij-J 

^0 



sin 



sin 



where 



and 



where 



c 

1 -f 2 cosfcuoT 

The last formula can be written in the more compact form 



pfe) = jW) cos — ^ + J™ sin — a i3 - 

C<J u 



d(»J') _ t(*J') ^ t(»J') 



r(»J') • U 
*2B Sin— «y, 

CJ 



AH) _ 
J la — 



■00 

00 



dTi 



J™ = /"°° dT sin ujT 

J — 00 



/-oo f sin 
JWJ = / dTcoswo^coso/N 

J— 00 




^sin(^f)' 



/ —00 



sin (^fj 



J25 = / c?T cos cj ^ sm S — 
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sin 


'2Rnu gin j 






sin (^f 


) 



cj T = + 27r/, v ? G(-7T,vr), Z = 0, ±1, ±2, 



(79) 



(80) 



(81) 

(82) 
(83) 

(84) 
(85) 
(86) 
(87) 

(88) 



we can transform the T-integral into the sum of the telescopic integrals according to the 
scheme: 



/oo 
dT 
-00 



/ — OO /* 



J 2=oo 
u <=— 00 

Using the fact that for the odd functions /(</?) and g(Z), the relations are valid 

1=00 



(89) 



f /(</?)^ = 0; J2 9(1) = 0, 



(90) 



i=— 00 



we can write 



UJq , J-TT I 



UJ UJ 1 | Sm 

OS (p COS 2717 

UJq UJq 



2Rnu> 



sin (f ) 



sin (f ) 



4 3) = o. 



For integrals with indices A, B we get: 



Aij) 
J 2A 



1 f T 

^2 COS if 

UJq i J—tt 



UJ w 1 I sm 

cos — ip cos 2irl — 
uj uj 



J 2B — u > 



2Rnu> 



sm 



sin (f ) 



So, the power spectral formula (80) is of the form: 



2 2 



where 



and 



pM = jM cos— ^ 

CJ 



-*2 — J 2A cos— a y . 



o; 



Using the Poisson theorem 



^ COS27T /= ^ UJq5(uJ — UJqI), 

l=-oo U k=-oo 



(91) 
(92) 

(93) 
(94) 

(95) 

(96) 
(97) 

(98) 



the definition of the Bessel functions J21 and their corresponding derivations and integrals 



1 

27 



dip cos ( z sin ^ ) cos lip = J21 (z) , 



— [ dipsva (zs\n sin(</?/2) coslip — — J^z), 
In J-tt \ 2 J 

1 riT sin (z sin |J /-z 

T" / — r V~7^ C0S ^= / J2i{x)dx, 
2tt J-tt sm{<p/2) Jo 



and using equations 



53 (-l) i+i cos — a - = 2(1 - cos 2ujt) = 4 sin 2 tut, 
ij=i w ° 

we get with the definition of the partial power spectrum Pi 

00 

P(u)=Y,6(u-luo)Pi, 
1=1 



(99) 
(100) 

(101) 
(102) 
(103) 



the following final form of the partial power spectrum generated by motion of two-charge 
system moving in the cyclotron: 



PAuj.t) = [4(sin ut) 2 }^^^ ( 2n 2 (3 2 J! 2l (2ln(3) - (1 - n 2 (3 2 ) dxJ 2l (x)) . (104) 

i\ri l v \ Jo J 

So we see that the spectrum generated by the system of electron and positron is formed 
in such a way that the original synchrotron spectrum generated by electron is modulated 
by function 4sin 2 (u;t). The derived formula involves also the synergic process composed 
from the synchrotron radiation and the Cerenkov radiation for electron velocity v > c/n 
in a medium. 

Our goal is to apply the last formula in situation where there is a vacuum. In this 
case we can put fj, — 1, n — 1 in the last formula and so we have 

Pi{u,t) = 4 sin 2 (ut) (2(3 2 4,(21(3) - (1 - (3 2 ) dxJ 2l (x)) . (105) 

7T V \ JO ) 

So, we see, that final formula describing the opposite motion of electron and positron 
in accelerator is of the form 

Pl,pair(u, t) = 4 Sin 2 (u)t) P^electron) (w) , (106) 

where P e i ec tron is the spectrum of radiation only of electron. For the same charges it is 
necessary to replace sine by cosine in the final formula. 

The result (106)is surprising because we naively expected that the total radiation of 
the opposite charges should be 

Pi(oj,t) = Pi 

(electron) y^i t) + Pl(positron) (w,t) . (107) 

So, we see that the resulting radiation can not be considered as generated by the 
isolated particles but by a synergic production of a system of particles and magnetic 
field. At the same time we cannot interpret the result as a result of interference of two 
sources because the distance between sources radically changes and so, the condition of 
an interference is not fulfilled. 

The classical electrodynamics formula (106) changes our naive opinion on the electro- 
dynamic processes in the magnetic field. From the last formula it follows that at time 
t = ixkjuj there is no radiation of the frequency uj. The spectrum oscillates with frequency 
uj. If the radiation were generated not in the synergic way, then the spectral formula would 
be composed from two parts corresponding to two isolated sources. 



5 The two center circular motions 

The situation which we have analyzed was the ideal situation where the angle of collision 
of positron and electron was equal to n. Now, the question arises what is the modification 
of a spectral formula when the collision angle between particles differs from tt. It can be 
easily seen that if the second particle follows the shifted circle trajectory, then the collision 
angle differs from ir. Let us suppose that the center of the circular trajectory of the second 
particle has coordinates (a,0). It can be easy to see from the geometry of the situation 
and from the plane geometry that the collision angle is tt — a, a ~ tana a/R where 
R is a radius of the first or second circle. The same result follows from the analytical 
geometry of the situation. 

While the equation of the first particle is the equation of the original trajectory, or 
this is eq. (70) 



xi(t) = x(t) = R(i cos(cj t) + j sin(cj t)), (108) 
the equation of a circle with a shifted center is as follows: 

x 2 (t) = x(i) = + cos(-cjot)) + j sm(-cj t)) = x(-t) + ia. (109) 

The absolute values of velocities of both particles are equal and the relation (72) is 
valid. Instead of equation (74) we have for radius vectors of particle trajectories: 

x^O-x^O^B^, (110) 

where Bn = An, B i2 = A 12 - ia, B 2 i = A 2 i + ia, B 22 = A 22 . 

In general, we can write the last information on coefficients By as follows: 

B i:j = A i:j + Eijia, (111) 

where e u = 0, e 12 = -1, £21 = 1, £22 = 0. 

For motion of particles along trajectories the absolute value of vector Ay ^> a during 
the most part of the trajectory. It means, we can determine Bij approximatively. After 
elementary operations, we get: 

|Bijl = ( A ij + 2|A ij |e ij acos^ J - + a 2 e\) xt2 , (112) 
where cos ipij can be expressed by the ^-component of vector and |A^-| as follows: 

cosy^ = ^f. (113) 

\ A ij\ 

After elementary trigonometric operations, we derive the following formula for (Aij) x : 

{A ij ) x = 2Rsm sm ~- ( 114 ) 

Then, using equation (114), we get with e = a/R 

\B i} \ = 2i?(sin 2 ^l±^l + ££i . sin ^ot+^or ^ c^r + ^e 2 _ ^ 

In order to perform the r-integration the substitution must be introduced. However, 
the substitution uqt + = ujqT does not work. So we define the substitution r = r(T) 
by the following transcendental equation (we neglect the term with e 2 ): 



. 2 UJ T + . U t + U T . UJ T 

sin + eeij sin sin 

o l J o o 



= sin— . (116) 

Or, after some trigonometrical modifications and using the approximative formula 
1 + x) 1/2 « 1 + x/2 for a; < 1 



[./.] 1/2 - sin ("° T 2 + ^) + |^sin ( 2 " 0T + 2 2 " 0t ~^) = sin^. (117) 

We se that for e = the substitution is c^ot + ctij = ^qT. The equation (117) is the 
transcendental equation and the exact solution is the function r = r(T). We are looking 
for the solution of equation (117) in the approximative form using the approximation 
sin a; ~ x. 

Then, instead of (117) we have: 



lu t + a 



e ( 2lo t + 2uj t — a tJ 



u T 



Using substitution 



uj t + = u T + u eA 
in eq. (118) we get, to the first order in e-term: 

Si 



2u 



l -^(uj T-2a ij + 2u;ot). 



Then, after some algebraic manipulation we get: 



c. 

uj t + aiij = uj T(1 - -Sij) - eeijUJot(-iy +J 



and 



ujt = uT(l - -Sij) (eeiji-iy^uot + otij) . 

2 Uq v ' 



For small time t, we can write approximately: 



and from eq. (122) 



cos(a> T + aiij) ~ cosu;o7 1 (l e^) 

2 



dr = dT(\ - -s i:j ) 



(118) 



(119) 



(120) 



121) 



(122) 



(123) 



(124) 



So, in case of the eccentric circles the formula (118) can be obtained from non- 
perturbative formula (80) only by transformation 



T 



T (! - a ij 



(££y(-l)* + W + aij) = aij, (125) 



excepting specific term involving sine functions. 
Then, instead of formula (80) we get: 



P(u,t) 



x cos(cjT j 

LO 



'4^2Rrf 
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(126) 



where T = T(l — §£ij). We see that only f and the a term are the new modification of 
the original formula (80). 

However, because e term in the sine functions is of very small influence on the behavior 
of the total function for finite time t, we can neglect it and write approximatively: 



P(u,t) = 



to e /i 
'4^2Rtf 



dT ]T {-iy +j 



x cos(cjT otij) 
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n- 



-v 2 cos(a;o7 1 ) 



sin 



2Rnu) 



sin 



( "oT 
\ 2 



sin (^) 



(127) 



So, we se that only difference with the original radiation formula is in variable otij. 
It means that instead of sum (102) we have the following sum: 



(-l) i+j cos— ocij = 2(1 - cos2cc>tcosac>t). (128) 
ij=i u ° 

It means that the one electron radiation formula is not modulated by [sinu;t] 2 but by 
the formula (128) and the final formula of for the power spectrum is as follows: 

Pi(u,t) = 2(1 - cos 2ut cos eut)Pi( e i ectr on ) (u) . (129) 
For e — * 0, we get the original formula (106). 

6 SUMMARY AND DISCUSSION 

We have derived, in the first part of the article, the total quantum loss of energy of the 
binary. The energy loss is caused by the emission of gravitons during the motion of the 
two binary bodies around each other under their gravitational interaction. The energy- 
loss formulas of the production of gravitons are derived here in the source theory. It is 
evident that the production of gravitons by the binary system is not homogenous and 
isotropical in space. So, the binary forms the "gravitational light house" where instead 
of the light photons of the electromagnetic pulsar are the gravitons. The detector of the 
gravitational waves evidently detects the gravitational pulses. 

This section is an extended and revised version of the older author's article (Pardy, 
1983a) and preprints (Pardy, 1994a,b), in which only the spectral formulas were derived. 
Here, in the first part of the article, we have derived the quantum energy-loss formulas 
for the linear gravitational field. Linear field corresponds to the weak field limit of the 
Einstein gravity. 

The power spectrum formulas involving radiative corrections are derived in the 
following part of this article, also in the framework of the source theory. The general 
relativity necessarily does not contain the method how to express the quantum effects 
together with the radiative corrections by the geometrical language. So, it cannot give 
the answer on the production of gravitons and on the graviton propagator with radiative 
corrections. This section therefore deals with the quantum energy loss caused by the 
production of gravitons and by the radiative corrections in the graviton propagator in 
case of the motion of a binary. 

We believe the situation in the gravity problems with radiative corrections is similar to 
the QED situation many years ago when the QED radiative corrections were theoretically 
predicted and then experimentally confirmed for instance in case of he Lamb shift, or, of 
the anomalous magnetic moment of electron. 

Astrophysics is, in a crucial position in proving the influence of radiative corrections 
on the dynamics in the cosmic space. We hope that the further astrophysical observations 
will confirm the quantum version of the energy loss of the binary with graviton propagator 
with radiative corrections. 

In the last part of this article on pulsars we have derived the power spectrum formula of 
the synchrotron radiation generated by the electron and positron moving at the opposite 
angular velocities in homogeneous magnetic field. It forms an analogue of the author 
article (Pardy, 1997) where only comoving electrons, or positrons was considered, and it 
forms the modified author preprints (Pardy, 2000a; 2001) and articles (Pardy, 200b; 2002,) 
where the power spectrum is calculated for two charges performing the retrograde motion 
in a magnetic field. The frequency of motion was the same because the diameter of the 
circle was considered the same for both charges. The retrograde motion with different 
diameters was not considered. 



It is surprising that the spectrum depends periodically on radiation frequency u> and 
time which means that the system composed from electron, positron and magnetic field 
behaves as a pulsating system. While such pulsar can be represented by a terrestrial 
experimental arrangement it is possible to consider also the cosmological existence in 
some modified conditions. 

To our knowledge, our result is not involved in the classical monographs on the 
electromagnetic theory and at the same time it was not still studied by the accelerator 
experts investigating the synchrotron radiation of bunches. This effect was not described 
in textbooks on classical electromagnetic field and on the synchrotron radiation. We hope 
that sooner or later this effect will be verified by the accelerator physicists. The radiative 
corrections obviously influence the synergistic spectrum of photons (Pardy, 1994c, d). 

The particle laboratories used instead of the single electron and positron the bunches 
with 10 10 electrons or positrons in one bunch of volume 300/im x 40/im x 0.01 m. So, 
in some approximation we can replace the charge of electron and positron by the charges 
Q and -Q of both bunches in order to get the realistic intensity of photons. Nevertheless 
the synergic character of the radiation of two bunches moving at the opposite direction 
in a magnetic field is conserved. 
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